Preprint typeset in JHEP style - HYPER VERSION 



On the supersymmetric vacua of the 
Veneziano-Wosiek model 



Matteo Beccaria 

Dipartimento di Fisica, Universita' di Lecce, Via Arnesano, 73100 Lecce 
INFN, Sezione di Lecce 



E-mail: matteo . beccaria@le . inf n . it 



ABSTRACT: We study the supersymmetric vacua of the Veneziano-Wosiek model in sec- 
tors with fermion number F = 2, 4 at finite 't Hooft coupling A. We prove that for F = 2 
there are two zero energy vacua for A > A c = 1 and none otherwise. We give the an- 
alytical expressions of both vacua. One of them was previously known, the second one 
is obtained by solving the cohomology of the supersymmetric charges. At F = 4 we 
compute the would-be supersymmetric vacua at high order in the the strong coupling ex- 
pansion and provide strong support to the conclusion that A = 1 is a critical point in this 
sector too. It separates a strong coupling phase with two symmetric vacua from a weak 
coupling phase with positive spectrum. 
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1. Introduction 

The study of quantum supersymmetric models for large number N of degrees of freedom 
has deep motivations in modern theoretical physics. Well known examples are matrix 
model formulations of M-theory [[|] and AdS/CFT duality between hf = 4 super Yang- 
Mills and type IIB superstring on AdS 5 x S 5 [|]. 

As illustrated in the recent review [||], direct techniques are currently available to 
analyze models in this class and a great deal of information can be obtained by combined 
analytical and numerical methods. These are based on an effective truncation of the state 
space with the minor drawback of introducing a controlled, and eventually irrelevant, 
supersymmetry breaking. Several examples at finite N are discussed in [^p. 

These methods can be extended to the most interesting limit N — > oo as explained in 
the beautiful series of papers [||, |6[ [7|, §J. Veneziano and Wosiek introduce a toy model of 
(non gauged) supersymmetric quantum mechanics at large N and show that the N — > oo 
limit can be described in terms of a planar Hamiltonian acting on single trace states. The 
dynamics in the planar limit is greatly simplified and non-trivial analytical and numerical 
results can be obtained (see also [||[l(]] for related developments). 

The Veneziano-Wosiek model is described in terms of N x N matrix fermion and 
boson creation/ annihilation operators with (non trivial) algebra 

[aij,a M ] = dudjk, {fij,fki) = s n?>jk- (1.1) 
Supersymmetry is generated by the nilpotent charges 

Q = Tr[/at(l + 5 at)], Qt = Tr[/ t a (1 + g a)), (1.2) 
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where g is a finite N coupling constant. The supersymmetric Hamiltonian is 

H = {Q,Ql}. (1.3) 



It commutes with Q, Qt, as well as with the additional operator C = [Q\Q] obeying 
C 2 = H 2 . The total fermion number F = Tr(/t /) is conserved. The total boson number 
B = Tr(at a) varies by AB = 0, ±1 under applications of H. 

In the large N limit, the Hamiltonian leaves invariant the subspace generated by sin- 
gle trace states of the form 

|m, . . . , n P ) = Tr[(at)"i f ■ ■ ■ (a*)"' /t] |0). (1.4) 

In this limit, the natural coupling turns out to be the 't Hooft combination A = g 2 N that 
will be kept fixed as N — > oo. 

The detailed analysis of the Veneziano-Wosiek model can be done in sectors with 
fixed F. At each F, one expects to find a certain number of supersymmetric vacua &f(A). 
The positive energy states are paired with supersymmetric partners in other sectors by 
the ladder action of the supersymmetry charges. The analysis of [||, 0, ||] investigates in 
details the F = 0, 1, 2, 3 cases at finite A and the general F sector at infinite A. The model 
turns out to be highly non trivial as we now summarize. 

The simplest sectors F = 0, 1 can be treated analytically. The number of supersym- 
metric vacua is 

6o(A) = a > i ' b ^=°- (l-^) 

The positive energy levels are almost evenly spaced and perfectly paired between the 
F = 0, 1 sectors by the action of the supersymmetry charges. At the critical point A = 1 
the spectrum collapses to zero. Across A = 1, there is a dynamical rearrangement of 
the SUSY multiplets. As a special feature, there is an exact strong/ weak duality holding 
separately in both sectors. 

The analysis of the next F = 2, 3 sectors is mainly numerical. At fixed A, the total 
number of boson excitations is truncated below a certain S max and the Hamiltonian is 
diagonalized. The spectrum is then extrapolated to B max — > oo limit. From the analysis, 
there are strong indications that 

HX) = { 2, A>t' ' hW = °> (L6) 

with a critical A c ~ 1. The positive energy levels are definitely not evenly spaced and they 
are only partially paired by the action of the supersymmetry charges. In other words, 
there are states with F = 3 which are not annihilated by Q^. They should be paired with 
states with F > 3 instead of F = 2. There is no sign of any weak/ strong coupling duality. 
As a general fact, the convergence of the numerical extrapolation f? max — > oo worsens as 
the critical point is approached. In principle, this can be a practical difficulty in obtaining 
accurate estimates for A c . 
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The sectors with F > 3 have been studied analytically by going to the extreme strong 
coupling limit A = oo where the boson number B is also conserved. This limit is partially 
solvable by mapping the Veneziano-Wosiek model to other models, i.e. a gas of g-bosons 
and, notably, the integrable XXZ spin 1/2 chain with anisotropy A = — ~. Many exact 
properties of this spin chain are known JTl|]. In particular, it is possible to predict the 
number of supersymmetric vacua. It reads 

Moo) = { 2,ifFG2N ' B = F±1 ' (1.7) 
1 0, otherwise. 

In the so-called magic case, F £ 2N and B = F ± 1, one can start from the two supersym- 
metric vacua at A = oo and write formal power series in 1/ y/X providing two zero energy 
states at finite A. For one of the two would-be vacua at F = 2 it can be shown that the 
state is normalizable for A > 1. For the other finite A vacuum at F = 2 and for the vacua at 
higher F the coefficients of the strong coupling series are known in a quite implicit form. 
The extension of the strong coupling phase is thus unknown. 

The supersymmetric vacua are states (p obeying the two equations 

Q<p = Qi<p = 0. (1.8) 

Due to Q 2 = (Q^) 2 = 0, they compute the cohomology of Q in the sector with F fermions. 
One can expect to take some advantage in determining zero energy states by solving the 
above pair of equations instead of solving directly the equation H tp = 0. 

In this paper, we follow this approach and extend the knowledge about bp(X) in two 
directions. First, at F = 2, we determine the all-order expression of the second supersym- 
metric vacua proving rigorously that A c = 1 in that case too. This completes the analysis 
of the F = 2 vacuum sector. Second, at F = 4 we compute at high order the strong 
coupling expansion of the vacua providing strong and accurate numerical support to the 
conclusion that again A c = 1. 

The plan of the paper is the following. In Sec. @, we give all the relevant formulae to 
work out the N — > oo planar limit of the Veneziano-Wosiek model. In Sec. @, we provide 
the exact analytical expressions of the F = 2 vacua. In Sec. (|J), we extend the analysis to 
the F = 4 sector. Sec. (@) is devoted to conclusions. 
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2. The N — > oo limit of the Veneziano-Wosiek model 



2.1 Hilbert space and norms 

In the N — > oo limit, the Hilbert space can be truncated to the H -invariant subspace gen- 
erated by single trace states of the form 

|n) = |n 1) ... 1 n F )=Tr[(atr/t...( a t)^/t]|0} 1 n t 6 N. (2.1) 

The basis states |n) obey 

|5n} = (-l) i7+1 |n), (2.2) 
where S is the left shift operator acting on sequences as 

5(ni, . . . ,n F ) = (n 2 , . . .,n P| ni). (2.3) 

The Hilbert space in the sector with F fermions is obtained by modding out the Zp action 
of S. The |n) states are not normalized. Their norm can be computed by applying the 
rules of planar calculus as explained in The result is 

|| |n) || 2 = d n N ni+ - +nF+F + subleading terms. (2.4) 

The multiplicity d n can be computed by the formula 

dn=E(- 1 )' (F+1) ^ n - (2-5) 

Of course, null states with d n = must be removed. With self-explanatory notation, the 
first cases are 

F = l:da = l, 

F = 2 : d aa = 0, else d ab = 1, 

F = 3 : d aaa = 3, else d abc = 1, 

F = 4 : d aa aa = 0, d a bab = 2 (fl ^ b), else d a bcd = 1, (2.6) 

F — 5 . d aaaaa — 5, else d a bcde — 1> 

F = 6 : daaaaaa = 0) d a 6aia6 = 3 (a ^ 6), d a bcabc = 0, else d a bcdef = 1- 

In the following we shall arbitrarily choose a representative in each 5-orbit and denote 
the resulting quotient Hilbert space as Hp. A simple choice amounts to lexicographically 
order n and shift it by S until the first element is mini n-i- I n particular, F = 2 basis states 
take the form |m, ng) with ni <ri2- 

2.2 Supersymmetry charges and Hamiltonian 

Let us split the supersymmetric charges in terms of operators with a definite variation of 
the boson number B. We write Q = Qi + Q2 with 

Qi= Tr[/at], (2.7) 
Q 2 = 5 Tr[/(at) 2 ], (2.8) 
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and similarly for the adjoint charges. We can go to the planar limit and remove all com- 
mon powers of N. For instance, the norms can be computed at leading order simply as 

II |n> || 2 = dn. (2.9) 

Then, by applying planar calculus, we obtain the following explicit formulas where the 
relevant coupling is indeed the 't Hooft combination A 

ni— 1 

Q\ \n u .. .,n F )= \k,ni - k - l,n 2 , . . .) - ^ \n lt k,n 2 - k - 1, n 3 , . . . ) + 

fc=0 fc=0 
"3-1 

\n 1 ,n 2 ,k,n 3 - k - 1, n 4 , . . . ) . (2.10) 

fc=0 

^ ni— 2 "2— 2 

-^=Q^|ni,...,n F ) = ]T | fc, m - A: - 2, n 2 , . . . ) - ^ |ni,fc,n 2 - A; - 2, n 3 , . . . } + 

"3-2 

X] 1^1,712, fc,n 3 - fc-2,7i4,...) , (2.11) 

fc=0 

and 

Qi |ni, . . . = |ni + n 2 + l,n 3 , . . . ) - \n u n 2 + n 3 + 1, n 4 , . . . ) + (2.12) 
\n 1 ,n 2 ,n 3 + n 4 + 1, . . . ) h (-l) F+1 |ni + n F + l,n 2 , . . . , nj?_i>. 

-^g 2 |ni,. . .,n F ) = |ni + n 2 + 2, n 3) . . . ) - |ni,n 2 + n 3 + 2,n 4 , . . . ) + (2.13) 

\ni,n 2 ,n 3 + n 4 + 2, . . . ) h (-l)^ 1 ]^ + n F + 2,n 2 , . . . ,n F -i). 

In principle, one could also write an explicit expression for the Hamiltonian. We shall not 
need it, apart from the F = 2 case. The expression is a bit involved and reads 

H\m,n 2 ) = [(m +n 2 + 2)(1 + A) - A(2 - S ni ,o + 26 ni+1>n2 )] \n u n 2 ) + 

+v / A [(m + 2) |m + 1, n 2 ) + (m + 1) \n x - 1, n 2 ) + (2.14) 
+(n 2 + 2) \m, n 2 + 1) + (n 2 + 1) |m, n 2 - 1)] + 

+2 A [(1 - <5ni + l,n 2 ) 1^1 + 1, n 2 - 1) + (1 - Jni.na + l) |™1 - 1, n 2 + 1)]. 

Notice that it involves A , A 1 / 2 , A 1 contributions. 
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3. Analytical supersymmetric vacua at F = 2 

From the explicit form of Q we obtain Q 7-L 2 = 0- This is in agreement with the SUSY 
multiplet structure of the F = 0, 1 sectors. As we recalled in the Introduction, all states in 
T-Li are in the image of Q^. Hence, given <p £ T-L 2 we have Q <p = Qt rj for some r] £ Hq. 
Hence, 

Q^Q<p = {Q^) 2 V = 0. ( 3. 1} 

This means that || Q ip || 2 = 0, or Q ip = 0. Therefore, the zero energy states in W2 are the 
solutions to the single equation 

Q 1 <P = 0. (3.2) 
To solve it, it is convenient to filter ip according to the boson number B. We write 

V = VBi (3-3) 

B>t 

where <ps has B boson excitations (there are no states with B = 1 in % 2 )- Replacing in 
Eq. Q, we find 

Q[ <p B + Ql <p B+1 = 0. (3.4) 
This equation determines ips+i in terms of ips as we now show. 



At B = 1, Eq. ( p^4[ ) has the unique solution 

pi = a|0,l). (3.5) 
At S = 2 we also obtain a unique solution 

^ 2 = -^|0,2). (3.6) 
The most general solution at B = 3 is 

<P3 = j |0,3}+/3 (-i|0,3) + 1 1,2)). (3.7) 

The existence of an additional arbitrary constant is consistent with the A = 00 analysis 
which also predicts that for B > 3 there should be no additional arbitrary constants. 
Taking ft = and a = 1 and iterating we easily obtain the explicit expression of the first 
supersymmetric vacuum 

This is normalizable for A > 1 and is the vacuum already found in 



-6- 



Taking a = and /3 = 1 we obtain the following expressions for the various terms of the 
second vacuum 



<Pi = f>2 = 0, 

-i|0,3) + |l,2) 



<P6 



2 

-L_(|0,4}-||1,3}), 

i(-^|0,4> + 5|M, + A| 2 ,3)), 

1/7 7 7 \ 



(3.9) 
(3.10) 

(3.11) 
(3.12) 
(3.13) 



and so on. Working out several additional levels, we are led to the following conjectured 
expressions 



r B - 1 



N B 



B-3 



, B 6 2N + 1, 



71 = 



(3.14) 



--1, 



with the coefficients 



CB,0 = ("I) 



B 7(B-1)(B-2)(B + 12) 



c B ,n = -210 (-1) 



2 (B + 2)(B + 3)(S + 4) 
s 5 - 2n 



(B + 2)(5 + 3)(B + 4) : 



B > 1, 
B > 2n. 



(3.15) 



We now prove that Eqs. (0.15 ) are indeed the unique solution of the basic equation 
Eq. fl3.2p . This is not as trivial as it could appear at first sight. Indeed, acting with 
produces states in H3 which are not necessarily in canonical order and the check is annoy- 
ing. An alternative proof exploits the more complicated equation H (p = which is more 
convenient in order to check the conjecture. Applying H to tp we obtain 

E TT^b- { ^ B + 2X 1 + A) - A (2 - <5„, - 2 S n+1 , B -n)] \n, B - n>+ 

\/A[(n + 2)|n + 1,S - n) + (S - n + 2)|n, B - n + 1) + (3.16) 

+(n + l)|n - 1, S - n) + (S - n + l)|n,S - n - 1)] + 

+2 A[(l - 5 n+x , B -n)\n + 1, B - n - 1) + (1 - <J n , B -n+i)|n - 1, 5 - n + 1)] } . 

This gives terms proportional to A or A leading to the two recursion equations 

{B + 5 n>0 - 2 5 n+ i,B-n) CB,n + (n + 1) c B -i,„-i + (S - n + 1) c B -i, n + (3.17) 

+2 (1 — <5 n ,B-n+l) CB,n-l +2(1— Jn+^B-n) CB,n+l = 0, 

(5 + 2) c Bin + (n + 2) cu+i.n+i + (5 - n + 2) c B+ i, n = 0. (3.18) 
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In the above two equations one has to set CB, n = if n does not obey < n < Nb- The 
boundary conditions that fix uniquely the solution to the recursion are 

c B , n = 0, £ = 1,2, (3.19) 
c 3 ,o = -\, (3.20) 
c 3 ,i = 1. (3-21) 

It is an easy check to verify that Eq. ( p.15 ) indeed solve the above recursion with the 
assigned boundary conditions. 

Thus, we have found the explicit expression for the second vacuum and it reads 

^ = E A ( B -3, /2(B + 2)( ( ^ > 3) ( B + 4) O - ^ - 2 '< B + 12 'l°' B > <M2 » 

N B 

-21Q^T{B -2ri)\n,B -n) 

n=l 

Evaluating the norm using the relation (holding MB £ N) 

N B i 

Y J {B-2n) 2 = -B(B-l)(B-2), (3.23) 

n=l 



we find 

II M (2) i|2 = « V ( g ~ 1) Cg ~ 2) (£ 3 + 19 B 2 + 36 B + 144) 

4 ^A A ' (B + 2)(B + 3) 2 (B + 4) 2 ^ ; 

The series converges for x = v < 1 and can also be resummed with the explicit result 

„ ll2 49 a; 4 + 20 z 3 - 1610 a; 2 - 8670 a; + 10260 
«W 2 = 1_ / 3 2 5) 

11 4 ai 6 (ai-l) v ' 

-73 ^ X ~ m l 1 + 7x) log(l - a;) - 44100^ U 2 (x), (3.26) 
a;' a;' 

where hi2(x) is the dilogarithm function. From this expression, we extract the singular be- 
havior in the x — > 1~ limit (the prefactor is an arbitrary normalization constant following 
our choice /3 = 1) 

ll*> (2) ll 2 = X7^ + ---- < 3 - 27 ) 

4 1 — x 

In conclusion, we have shown that the cohomology of Q at F = 2 has two solutions both 
normalizable for A > A c = 1 only i.e. 

M*) = { 2' a > i' ' < 3 ' 28 > 
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4. Results at F = 4 



4.1 A view to the spectrum 

The sector F = 4 is much more complicated. As a first step, we have diagonalized H up 
to B max = 24 to have a feeling about would-be zero modes in the S max — > oo limit. The 
smallest 6 levels are shown in Fig. @ for the values A = 0.5, 1, 1.2, and 4. For A = 0.5 and 
A = 4 it seems quite clear that there are respectively and 2 supersymmetric vacua. At 
A = 1 it is plausible that all levels are converging to zero in agreement with the reasonable 
conjecture that the critical point is again A c = 1. However, the estimate of A c is difficult at 
these values of S max as illustrated by the inset at A = 1.2. Here a clean stabilization as for 
A = 4 would require quite a larger i? max . If we do not want to push further the numerical 
diagonalization, it seems mandatory to find an alternative determination of the critical 
point. 

Are the methods exploited at F = 2 applicable ? In the next Section we shall address 
this question discussing some difficulties and their (numerical) resolution. 

4.2 Strong coupling expansion of Q cohomology 

At F = 4 we find zero energy states by imposing the full set of cohomological equations 

Q<p = Ql<p = 0. (4.1) 

The general solution is quite complicated compared to the case F = 2 and the solution 
does not organize well in powers of A. The reason is that ip produces a series in de- 
scending powers of V\ as at F = 2. However, the equation Q <p = has the opposite 
behavior. 

We can bypass this problem recalling that, after all, we are interested in the determi- 
nation of the convergence radius of the strong coupling expansion. Thus, we try to solve 
Eqs. ([P]) by making from start the Ansatz 

oo oo 

<P=T, A "" /2 E <Pn,B, (4.2) 
n=0 B=l 

where, as indicated, (p nt B is a state with boson number B. Removing the \/A factors in Q 2 
and Ql we have to solve the equations 

Q2 <fn,B + Ql <Pn-l,B+l = 0, (4.3) 
Q\ <Pn,B + Ql <Pn-l,B-l = 0. 

It is easy to check that these equations are compatible and admit a unique solution for 
the n-th order in terms of (p computed at (n — l)-th order. This is true with the exception 
of those values of B where the operator Q 2 has non empty cohomology. However, the 
cohomology of Q 2 is given by the zero energy states of the A = oo Veneziano-Wosiek 
model which is known. It contains a state at each B = F ± 1, here B = 3,5. The explicit 
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zero modes (with an arbitrary normalization) are 



(773) = |0, 0, 0, 3} - 3 |0, 0, 1, 2} + 3 |0, 0, 2, 1} + 4 |0, 1, 0, 2} - 7 |0, 1, 1, 1), (4.4) 
|77s> = -10,0,0,5) +410,0,1,4) -7|0,0,2,3) +7|0,0,3,2) -4|0,0,4,1) + (4.5) 
-6 |0, 1,0,4) + 1710,1, 1,3) - 21 10,1,2,2) + 14 |0, 1,3, 1) - 12 |0, 2, 0,3) + 
+25 |0, 2, 1, 2) - 21 10, 2, 2, 1) + 17 |0, 3, 1, 1) - 42 |1, 1, 1, 2). 

Let us discuss in some details the solution which reduces at A = 00 to 1 773). The other case 
is completely similar. We start with 

<Po = J2 <P°,b = 1^3) • (4-6) 

B 

Then, we solve at each B Eq. ( |Q| ) with n = 1. Of course, there is a maximum B beyond 
which we do not have non vanishing solutions for p 1B . The procedure is iterated. The 
solution of Eq. Q4.3P is always unique with the exception of the cases B = 3, 5 where we 
can add to <p n ,B an arbitrary constant a times \t]b)- The general solution can always be 
put in the form 

<Pn,z = <$° m + a\0,l,l,l), (4-7) 
Vn, 5 = ^° m + a|l,l,l,2). (4.8) 

In other words, the inhomogeneous piece of the solution does not have contributions from 
the states |0, 1, 1, 1) and |1, 1, 1, 2) which come totally from the zero modes. We arbitrarily 
set a = to fix the zero mode contributions. Other choices are possible, but do not change 
the convergence properties of the strong coupling expansion. 

The explicit expression of ip nt B is quite complicated and unfortunately we did not succeed 
in finding a closed formula. However one can try to estimate the convergence radius from 
a study of the strong coupling series. To this aim, we have evaluated the norm of the 
would-be vacuum by working out the terms of 



00 

<P II 2 = E (4-9) 



71 = 



up to n = 24. After normalization, the first terms read 



2 _ 33163 1 3606544643 1 6669989903943227 1 

^" ~~ + 13272 A + 777060669 X + 891227976610818 A^ + 
2155907292859955213802297145858 1 

+ 195526455552229171879002565071 A 1 + 
5974912975520703560602997490582425731877554327739 1 

+ 393252674274544022631209524089738568938469988544 A^ + 
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The convergence radius can be estimated by the ratio test or, better, by means of improved 
recurrent estimators in the spirit of [12]. In other words, we compute the sequences 



d(0) 
tin 


a 


+ 1 

i 

n 


■Tin 


= n 

Q-n 




1 

~ 2 


n 2 - 


p(3) 
tin 


1 

~ 6 


o a 
n d - 




^ 5^ " 



(n-1) 



3(n- l) 2 



0-n-\ 



a n -2 



+ 3(n-2) 2 



&n-2 



(n-3) 



3 a„_2 



j(3) 



(4.10) 
(4.11) 

(4.12) 

(4.13) 

(4.14) 



Fig. @ show the results obtained with rW, R( 2 \ and i? (3) . As one can see, the 
ratio test (sequence R^) is poorly useful in determining A c . Instead, the higher order 
estimators converge more and more quickly to a A c that can be estimated to be 



A c = 1.000(1). 



(4.15) 



This result shows that there is a supersymmetric vacuum extending up to A c ~ 1 at F = 4. 
The same procedure can be started from the strong coupling vacuum at B = 5, repeating 
the construction and removing the component along the first vacuum in order to enforce 
orthogonality. The numerics is less clean, but fully consistent with the above estimate. 
Thus, we have provided strong support to the conclusion that 



64(A) 



0, A < 1, 
2, A > 1 



(4.16) 



It is clear that the methods described in this section can be extended to larger F with no 
additional difficulties. 



5. Conclusions 

The Veneziano-Wosiek model is a surprisingly rich toy model where quantum supersym- 
metry at large N can be investigated. As is usual in super symmetry, a lot of information 
is already contained in the most basic question, the dimension of the vacuum sector, the 
integer number bp(X). In this paper, we have extended the known results for F = 0, 1, 2, 3 
providing new analytical results at F = 2 and F = 4. Our results support the conjec- 
ture that the two strong coupling supersymmetric vacua existing for even F > 2 can be 
analytically continued up to the critical value A c = 1 in all fermion sectors. 

Most interestingly, the Veneziano-Wosiek model is known to have some intriguing 
connection with combinatorial problems as discussed in [gj. This fact is well established in 
the extreme strong coupling limit. The mapping to the XXZ spin chain permits to extend 
to the Veneziano-Wosiek model several number-theoretical facts [ |TT| ] recently exploited 
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in the context of Alternating sign matrix conjectures p4fl. Similar relations between su- 
persymmetric models and combinatorics are actually not new as discussed in the SUSY 
algebra non-linear realizations discussed in [ fi3| ] and also related to the XXZ chain at the 
peculiar anisotropy A = — |. 

What is somewhat surprising is the fact that hidden combinatorial facts could be at 
work even at finite coupling. The search for supersymmetric vacua in the F = 2 sector 
described in this paper has been achieved due to the ability of guessing the solution of a 
complicated recursion problem. As soon as a guess is proposed, it can be checked with 
minor effort. However the guess itself was not trivial. Actually we could find it by search- 
ing within suitable classes of rational sequences arising precisely in typical combinatorial 
problems 



Acknowledgments 

We thank J. Wosiek for communications about his numerical results on the model and G. 
F. De Angelis for conversations. 



References 

[1] T. Banks, W. Fischler, S. H. Shenker and L. Susskind, M theory as a matrix model: A conjecture, 
Phys. Rev. D 55, 5112 (1997) [arXiv:hep-th/9610043]. 

[2] J. A. Minahan, A Brief Introduction To The Bethe Ansatz In M = 4 Super-Yang-Mills, J. Phys. A 
39, 12657 (2006). 

[3] J. Wosiek, Solving some gauge systems at infinite N, Lectures given at Cracow School of 
Theoretical Physics: 46th Course 2006, Zakopane, Poland, 27 May - 6 Jun 2006. 
[arXiv:hep-th/0610172]. 

[4] P. Bialas and J. Wosiek, Lattice study of the simplified model ofM-theory, Talk given at 36th 
Rencontres de Moriond on QCD and Hadronic Interactions, Les Arcs, France, 17-24 Mar 
2001, [arXiv:hep-lat/0105031]. 

J. Wosiek, Spectra of supersymmetric Yang-Mills quantum mechanics, Nucl. Phys. B 644, 85 
(2002) [arXiv:hep-th/0203116]. 

J. Wosiek, Supersymmetric Yang-Mills quantum mechanics, Invited talk at NATO Advanced 
Research Workshop on Confinement, Topology, and other Nonperturbative Aspects of 
QCD, Stara Lesna, Slovakia, 21-27 Jan 2002, arXiv:hep-th/0204243. 

M. Campostrini and J. Wosiek, Exact Witten index in D = 2 supersymmetric Yang-Mills 
quantum mechanics, Phys. Lett. B 550, 121 (2002) [arXiv:hep-th/0209140]. 
M. Trzetrzelewski and J. Wosiek, Quantum systems in a cut Fock space, Acta Phys. Polon. B 
35, 1615 (2004) [arXiv:hep-th/0308007]. 

J. Wosiek, Recent progress in supersymmetric Yang-Mills quantum mechanics in various 
dimensions, Acta Phys. Polon. B 34, 5103 (2003) [arXiv:hep-th/0309174]. 

M. Campostrini and J. Wosiek, High precision study of the structure ofD = 4 supersymmetric 
Yang-Mills quantum mechanics, Nucl. Phys. B 703, 454 (2004) [arXiv:hep-th/0407021]. 

J. Wosiek, Supersymmetric Yang-Mills quantum mechanics in various dimensions, Int. J. Mod. 
Phys. A 20, 4484 (2005) [arXiv:hep-th/0410066]. 



-12- 



J. Wosiek, On the SO(9) structure of supersymmetric Yang-Mills quantum mechanics, Phys. Lett. 
B 619, 171 (2005) [arXiv:hep-th/0503236]. 

J. Wosiek, Vacua of supersymmetric Yang-Mills quantum mechanics, Contributed to 11th 
International Conference on Elastic and Diffractive Scattering: Towards High Energy 
Frontiers: The 20th Anniversary of the Blois Workshops, Chateau de Blois, Blois, France, 
15-20 May 2005. [arXiv:hep-th/0510025]. 

[5] G. Veneziano and J. Wosiek, Planar quantum mechanics: An intriguing supersymmetric 
example, JHEP 0601, 156 (2006) [arXiv:hep-th/0512301]. 

[6] G. Veneziano and J. Wosiek, Large N, supersymmetry ... and QCD, Sense of Beauty in Physics 
- A volume in honour of Adriano Di Giacomo, edited by M. D'Elia, K. Konishi, E. 
Meggiolaro and P. Rossi (Ed. PLUS, Pisa University Press, 2006), [arXiv:hep-th/0603045]. 

[7] G. Veneziano and J. Wosiek, A supersymmetric matrix model. II: Exploring 
higher-fermion-number sectors, JHEP 0610, 033 (2006) [arXiv:hep-th/0607198]. 

[8] G. Veneziano and J. Wosiek, A supersymmetric matrix model. Ill: Hidden SUSY in statistical 
systems, JHEP 0611, 030 (2006) [arXiv:hep-th/0609210]. 

[9] E. Onofri, G. Veneziano, J. Wosiek, Supersymmetry and Combinatorics, 
[arXiv:math-ph/0603082] . 

[10] R. De Pietri, S. Mori and E. Onofri, The planar spectrum in U(N)-invariant quantum mechanics 
by Fock space methods. I: The bosonic case, [arXiv:hep-th/0610045]. 

M. Bonini, G. M. Cicuta and E. Onofri, Fock space methods and large N, 
[arXiv:hep-th/0701076]. 

[11] F. C. Alcaraz and V. Rittenberg, Supersymmetry on Jacobstahl lattices, J. Phys. A 38, L809 
(2005) [arXiv:cond-mat/0510272]. 

J. de Gier, A. Nichols, P. Pyatov and V. Rittenberg, Magic in the spectra of the XXZ quantum 
chain with boundaries at Delta = and Delta = -1/2, Nucl. Phys. B 729, 387 (2005) 
[arXiv:hep-th/0505062]. 

A. V. Razumov, Yu. G. Stroganov, O(l) loop model with different boundary conditions and 
symmetry classes of alternating-sign matrices, Theor.Math.Phys. 142 (2005) 237-243; 
Teor.Mat.Fiz. 142 (2005) 284-292. 

A. V. Razumov, Yu. G. Stroganov, Combinatorial nature of ground state vector ofO(l) loop 
model, Theor.Math.Phys. 138 (2004) 333-337; Teor.Mat.Fiz. 138 (2004) 395-400. 

A. V. Razumov, Yu. G. Stroganov, Spin chains and combinatorics: twisted boundary conditions, 
J.Phys. A 34 (2001) 5335-5340. 

A. V. Razumov, Yu. G. Stroganov, Spin chains and combinatorics, J.Phys. A 34 (2001) 3185. 
Yu. Stroganov, The Importance of being Odd, J.Phys. A 34 (2001) L179-L186. 

[12] Y. F. Chang and G. Corliss, Ratio-Like and Recurrence Relation Tests for Convergence of Series, 
IMA Journal of Applied Mathematics 1980 25(4):349-359;doi:10.1093/imamat/25.4.349. 

Y. F. Chang and G. Corliss, Solving Ordinary Differential Equations Using Taylor Series, ACM 
Transactions on Mathematical Software (TOMS) archive, Volume 8 , Issue 2 (1982), 114 - 144 
ISSN:0098-3500. 



-13- 



[13] P. Fendley and K. Schoutens, Exact results for strongly-correlated fermions in 2+1 dimensions, 
Phys. Rev. Lett. 95, 046403 (2005) [arXiv:cond-mat/0504595]. 

P. Fendley K. Schoutens and B. Nienhuis, Lattice fermion models with super symmetry, J. Phys. 
A 36, 12399 (2003) [arXiv:cond-mat/0307338]. 

Xiao Yang, Paul Fendley, Non-local space-time supersymmetry on the lattice, J. Phys. A 37 
(2004) 8937. 

M. Beccaria and G. F. De Angelis, Exact Ground State and Finite Size Scaling in a 
Supersymmetric Lattice Model, Phys. Rev. Lett. 94, 100401 (2005) [arXiv:cond-mat/ 0407752]. 

[14] P. Di Francesco, A refined Razumov-Stroganov conjecture II, J. Stat. Mech. 0411, P004 (2004) 
[arXiv:cond-mat/0409576]. 

A. V. Razumov, Yu. G. Stroganov, Bethe roots and refined enumeration of alternating-sign 
matrices, J. Stat. Mech. (2006) P07004. 

A. V. Razumov, Yu. G. Stroganov, Enumeration of quarter-turn symmetric alternating-sign 
matrices of odd order, [arXiv:math-ph/ 0507003]. 

A. V. Razumov, Yu. G. Stroganov, Enumerations of half-turn symmetric alternating-sign matrices 
of odd order, Theor.Math.Phys. 141 (2004) 1609-1630; Teor.Mat.Fiz. 141 (2004) 323-347. 

[15] Christian Krattenthaler, Advanced Determinant Calculus, Sminaire Lotharingien de 
Combinatoire, B42q (1999), http : / /www .mat . univie . ac . at / sic/ 



-14- 



F = 4 




o.oi 



10 12 14 16 18 20 22 24 




10 12 14 16 18 20 22 24 



Figure 1: Smallest 6 energy levels at F = 4 and various couplings A as functions of the upper limit 
on the boson number B < B^^. 
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